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We determine the asymptotic level spacing distribution for the Laguerre En- 
semble in a single scaled interval, (0, s), containing no levels, Ep(0,s), via Dyson's 
Coulomb Fluid approach. For the a = Unitary-Laguerre Ensemble, we recover the 
exact spacing distribution found by both Edelman and Forrester, while for a / 0, 
the leading terms of £^(0, s), found by Tracy and Widom, are reproduced without 
the use of the Bessel kernel and the associated Painleve transcendent. In the same 
approximation, the next leading term, due to a "finite temperature" perturbation 
(/3 + 2), is found. 



Introduction 

The probability that there are no levels in a scaled interval (—t,t) [where t is mea- 
sured with respect to the averaged spacing], E(0,t), in a long stack of energy levels 
of heavy nuclei was given by a conjecture of Wigner|l| and was well supported by exper- 
imental data|2| for systems with time reversal symmetry. In a series of seminal papers, 
Dyson introduced a new class of random matrix ensembles and determined in the con- 
tinuum approximation (expected to be valid when number of levels, N, is very large) 

that In #0(0,*) \f3t 2 - (1 - (3/2)irt, using the methods of classical electrostatics, 

potential theory and thermodynamics, for ensembles with unitary (/3 = 2), orthogonal 
{(3 = 1) and symplectic (fl = 4) symmetries. A term of O(lnt) and a constant missed in 
the continuum approximation was later discovered by Widom, des Cloizeaux and Mehta, 
and Dyson § §, § 0, §. 

Recently, in a remarkable paper, Tracy and Widom|| showed, in the single interval case, 
that the logarithmic derivative of the Fredholm determinant of the Bessel Kernel — which 
arises in the scaling limit [with respect to the hard edge density, see[|| [H|] of the unitary 



Laguerre ensemble — satisfies a Painleve V equation, from which the asymptotic level 
spacing distribution can be computed exactly, amongst other quantities of interest from 



the random matrix point of view [1 'I 



In this letter we shall employ the continuum approximation of Dyson to calculate the 
level spacing distribution for the Laguerre ensemble [with a / 0, see||], with (3 = 2. 
For P 7^ 2, the spacing distribution can be found by a perturbative calculation due to 
Dyson§. 

From the Brownian motion model, in the "hydrodynamical" approximation, Dyson E 



derived an equation satisfied by the non-equilibrium level density a(x, r); 

d d ( 

s^^sHsJ' T>( > (1) 

with 

r) = - y dya(y, r) In |x - y\ + f ^ - ^ In r)] , (2) 

where the fictitious time r pulls the levels towards the observed level density, r(x), gen- 
erated by the imposed potential u(x)/(3 = J dyr(x) In \x — y\, (as r — > oo), that holds the 
Coulomb Fluid together. The stationary solution, reached as r — > oo, for the level density 
statisfies the Hiikel-like self-consistent equation 

u(x) — {3 J dya{y) In \x — y\ + ^1 — ^\ In [<j{x)\ = A = constant, (3) 

with effective temperatute T e oc (1 — f ). This, in turn, maybe derived from the following 
variational principle: min CTj ^ F[a, //], with 

F[a, n]= J dx^(x)cr(x) - fi(J dxa(x) - N^j , (4) 

and A = —\i + (1 — (3/2), where \x is the chemical potential. Therefore the free energy, F, 
at equilibrium, with exactly N levels contained in an interval I is 



F[I] = ^AN + i J^dx u(x)a(x) - ~ (l - J^dxa(x)\n[a(x)}, 



(5) 



subject to jjdxa{x) = N. For the Laguerre Ensemble, u(x) = x — alnx, x G (0, oo), but 
as we are using the continuum approximation, an upper band edge, b G (0, oo) must be 
imposed on the level density to produce a finite number of levels. In the large N limit, 



the level density is a(x) = -^y x x , % G (0,4iV) for the unitary case[13, 14]. This is 
distinct from the Wigner semi-circle law, where u(x) = x 2 , x G (— oo, +oo). 



Asymptotic Spacing 

The probability distribution that an interval (0, a) contains no levels is, by definition, the 
ratio of the partition function where all N levels reside in the complement of (0, a) i.e. 
(a, b) [where b is the upper band edge and < a < b] to that for which all N levels reside 
in the full interval i.e. (0, b), and is|§, 

ln^(0,a) = -[F(a,6)-F(0,6)]. (6) 

In the continuum approximation, F is given by Eq.(H), where cr(x) solves Eq. (RH [T5|. 



From thermodynamic considerations, since F(a, b) is the Free Energy in a "constricted" 
region, we must have F(a,b) > F(0,b). We now proceed to solve Eq. (||) for (3 = 2, 
in the interval (a,b). To simplify the mathematics, Eq.(|||) is converted into a singular 
intergral equation by taking a derivative with respect to x, and can be solved by a standard 



method 16 1. The solution of this equation is 



1 lb — x r b dy y — af a\ 1 b — x f a fa\ , . 

a[X) = ^x~^J a J^ib^y y 1 ~y) = ~^x~^~aV~xib) ' (?) 

where x G (a, b) and to maintain positivity we demand that \fab > a. A straightforward 
calculation supplies the normalization condition: 



b — a a 



1- 



We have deliberately left f3 without setting it equal to 2 in Eqs. (0) and @. It is clear from 
the structure of Eq.(^), that N, the total number of levels, is almost exhausted by the first 
term, but as to be seen later, the second term cannot be discarded. As we are required to 
compare F(a, b) with F(Q, b), where iV is very large, to facilitate the computation we shall 
evaluate instead, F(a,b) in the limits a < 6, rf- 1 and r^- ~ 0(1), thus by-passing 



an independent computation of F(0,b). To evaluate F(a,b) requires the determination 



of A and the "interaction energy" \ J a dxa(x)u(x). In the limits stated we shall extract 



the very large terms, which are functions of iV only and finite terms that are functions 
of Na only; any remaining terms are therefore negligable in the large N limit. Note that 
the very large terms are then subtracted according to the definition of Ep(0,a), leaving 
behind terms only those terms in Na. 

For the constant A, we send s — > b~ in Eq.([|) [with (3 = 2], which gives 



a b — a . 
A = a In 



4e 



a In b 



a 



In 6- a + -,/-/. . 

7T V \b 



(9) 



where 



Hx) 



V * i + a; 



J "5-") F ( 1 -5-5 + A -i 



-2tt(1 - In 2) + vrv^r H , x < 1. 



£=3/2 



(10) 



As N — > 00, but with iVa finite, we find that 



-NA N 2 ln(JV/e) - - ln(16iV/e 4 ) + ^Na - ^VNa 



'ir 



where to capture all the terms in Na it is essential to use Eq.(|) for N. For the contribution 
due to the interaction, we have (for (3 = 2) 



2 1 <*<*>«<«> = 5 



1 (6 — a) a (6 — a) 
+ 2~ a_ 2/3 2 2(3 



a [b — a) 1 a 
h 



7T 2(3 \b-aj 2(3 



a 



— ln(6 — a) 



I" 1 



a a 
+ —\It9 



2/3tt V 6 " \6- a 



12) 



where 



r 1 /l — £ 7r 

/i(x) := / dt 4/ ln(t + 2) ln(4e) - tt a; + 277^ + ■ 

Jo V t 2 



s«l (13) 



and 



r 1 dt 1-t, , 
9{x) := / — — \— — \n(t + x) 
Jo t + x V t 



ix d 
~2dp 



, -''/■'( p, i; 2; -1/z 



~ — p ln(4x) + 7rlnxH , sc < 1 (14) 

p=i v # 



With the same considerations, 

1 i-b i / a 2 1 3 / a 2 

-/ dta(t) ~ -iV 2 -aiVln(2/e 1/2 ) ln4JV + -Na- -av / iV^+— ln(4iVa). (15) 

2 Ja 2 2 2 2 4 

Pooling together Eqs. ( |TTj) and fll5|) we find by subtracting the very large terms 



a 2 



which gives, 



F(a, b) - F(a < b, b) ~ iVa - 2aViVa + — ln(4iVa), (16) 



3»(0, -to , (17) 



upon scaling with respect to the hard edge density, i.e. with the replacement Na —>■ |. 
Before we proceed to give the result for the (3^2 case, we should like to mention 
that by repeating the calculation in the far simpler situation where a = 0, Eq.(^) can 
be solved trivially, which gives b — a = 4N [for (3 = 2} and the change in the free 
energy is equal to Na or |. This is the exact result of Edelman|T7[] and of Forrester [IS|; 



E 2 (0,s) = e~ s / 4 , a = 0, s G (0,oo). For a ^ 0, (3 = 2, Tracy and Widom[g, found 
through an asymptotic expansion of a Painleve V equation, higher order terms, in l/\/s, 
1/s etc., and made a conjecture concerning the term independent of s. It was observed 
by Tracy and Widom that[[| with a = =f|, the Bessel kernel reduces to the kernels which 
arise by scaling into the bulk of the spectrum of the Gaussian orthogonal and symplectic 
ensembles, respectively, provided one makes the replacement s —>■ TT 2 t 2 . It is amusing 

6 



to see that by approching the GOE and GUE through the "back door", i.e. using the 

mapping of Tracy and Widom||; the Coulomb Fluid approach, when applied to the 

Laguerre ensemble, supplies rather precise information |19|| . 

Correction to the Free Energy When (3^2 

This is simply found by adding the finite temperature contribution 



7 — - \) = [ dxa(x) \n\a(x)] ~ — In 
(l-f) J* 



2tt 



a a 

N ln4iV+-ln4jVa, a<6-a, (11 

2 4 



to the f3 = 2 Free Energy. 

Collecting the appropriate terms together, we find that 

- In E a (0, s) ~ — - — + — In s + i^l - | J - In 8 , (19) 

which clearly reduces to Eq.(^) when (3 = 2. 
Conclusion 

We should like to mention that by treating the constraint, i.e. Eq.(§) more accurately, it 
should be possible to produce the higher order correction terms l/y/s, Ijs ... mentioned 
previously. However, the computation would become exceedingly complicated and clearly 
the methods of Tracy and WidomQ have much to be desired. 

The above calculations suggest upon comparison with exact results that the Coulomb 
Fluid approach is quite robust and may shed light on the level spacing distribution of the 
q— Laguerre ensemble, which arises in the context of electronic transport in disordered 
systems [0, with the potential 

oo 

M (x,g) = $>[l + (l-g)xgT, q G (0, 1) (20) 

n=0 

which reduces to the ordinary Laguerre potential as q — > 1~; u(x, 1~) = x. We leave as a 
future project the determination of the level spacing distribution with u(x; q). 
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